Intro

Issue

TROLL model currently compute leaf lifespan with Reich’s allometry (Reich et al. 1991). But we have shown that Reich’s allometry is underestimating leaf lifespan for low LMA species. Moreover simulations estimated unrealistically low aboveground biomass for low LMA species. We assumed Reich’s allometry underestimation of leaf lifespan for low LMA species being the source of unrealistically low aboveground biomass inside TROLL simulations. We decided to find a better allometry with Wright et al. (2004) GLOPNET dataset.

We tested different models starting from complet model Mcomp: \[ {LL_s}_j \sim \mathcal{logN}(log({\mu_s}_j),\,\sigma)\,, ~~s=1,...,S_{=4}~, ~~j=1,...,n_s\]

\[{\mu_s}_j = {\beta_0}*e^{{\beta_1}_s*{LMA_s}_j^{{\beta_3}_s} - {\beta_2}_s*{Nmass_s}_j^{{\beta_4}_s}}\] \[{\beta_i}_s \sim \mathcal{N}({\beta_i},\,\sigma_i)\,^I\] \[(\beta_i, \sigma, \sigma_i) \sim \mathcal{\Gamma}(0.001,\,0.001)\,^{2I+1}\]

LL graph

Figure 1: Leaf mass per area (LMA), leaf nitrogen content (Nmass) and leaflifespan (LL). Leaf mass per area (LMA in \(g.m^{-2}\)), leaf nitrogen content (Nmass, in \(mg.g^-1\)) and leaf lifespan (LL in \(months\)) are taken in GLOPNET dataset from Wright et al. (2004).

M1

Model

\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA + {\beta_2}_s*N,\sigma)\,\] Maximum likekihood of 11.0675641 and \(R^2\) of 4.011824110^{30}

Convergence

M2

Model

\[ LL \sim \mathcal{logN}(\beta_0 + LMA^{{\beta_3}_s} + N^{{\beta_4}_s},\sigma)\,\] Maximum likekihood of 8.7566631 and \(R^2\) of 1.612299710^{32}

Convergence

M3

Model

\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA^{{\beta_3}_s} + {\beta_2}_s*N^{{\beta_4}_s},\sigma)\,\] Maximum likekihood of 9.3754692 and \(R^2\) of 1.017729710^{31}

Convergence

M4

Model

\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA + N,\sigma)\,\] Maximum likekihood of 8.9827971 and \(R^2\) of 3.283734610^{24}

Convergence

M5

Model

\[ LL \sim \mathcal{logN}(\beta_0 + LMA^{{\beta_3}_s} + N^{{\beta_4}_s},\sigma)\,\] Maximum likekihood of 9.1974485 and \(R^2\) of 3.554504210^{27}

Convergence

M6

Model

\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA^{{\beta_3}_s} + N^{{\beta_4}_s},\sigma)\,\] Maximum likekihood of 9.257722 and \(R^2\) of 6.074544210^{27}

Convergence

M7

Model

\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA,\sigma)\,\] Maximum likekihood of 8.9051831 and \(R^2\) of 5.398465610^{22}

Convergence

M8

Model

\[ LL \sim \mathcal{logN}(\beta_0 + LMA^{{\beta_3}_s},\sigma)\,\] Maximum likekihood of 9.1917157 and \(R^2\) of 6.528413310^{24}

Convergence

M9

Model

\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA^{{\beta_3}_s},\sigma)\,\] Maximum likekihood of 8.931063 and \(R^2\) of 6.279900410^{29}

Convergence

Results

Results

M1 M2 M3 M4 M5 M6 M7 M8 M9
ML 1.106800e+01 8.7570e+00 9.37500e+00 8.983000e+00 9.197000e+00 9.258000e+00 8.905000e+00 9.192000e+00 8.9310e+00
R2 4.011824e+30 1.6123e+32 1.01773e+31 3.283735e+24 3.554504e+27 6.074544e+27 5.398466e+22 6.528413e+24 6.2799e+29

Figure 3: Model predictions.

References

Reich, P.B., Uhl, C., Walters, M.B. & Ellsworth, D.S. (1991). Leaf lifespan as a determinant of leaf structure and function among 23 amazonian tree species. Oecologia, 86, 16–24.

Wright, I.J., Reich, P.B., Westoby, M., Ackerly, D.D., Baruch, Z., Bongers, F., Cavender-Bares, J., Chapin, T., Cornelissen, J.H.C., Diemer, M. & Others. (2004). The worldwide leaf economics spectrum. Nature, 428, 821–827.